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Color screening in a constituent quark model of hadronic matter 
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The effect of color screening on the formation of a heavy quark-antiquark (QQ) bound state — 
such as the J/ip meson — is studied using a constituent-quark model. The response of the nuclear 
medium to the addition of two color charges is simulated directly in terms of its quark constituents 
via a string-flip potential that allows for quark confinement within hadrons yet enables the hadrons 
to separate without generating unphysical long-range forces. Medium modifications to the proper- 
ties of the heavy meson, such as its energy and its mean-square radius, are extracted by solving 
Schrodinger's equation for the QQ pair in the presence of a (screened) density-dependent poten- 
tial. The density dependence of the heavy-quark potential is in qualitative agreement with earlier 
studies of its temperature dependence extracted from lattice calculations at finite temperature. In 
the present model it is confirmed that abrupt changes in the properties of the J/V'-meson in the 
hadronic medium (plasma), correlate strongly with the deconfining phase transition. 

PACS numbers: 24.85.+p, 24.10.Lx, 25.75.Nq, 12.39.Jh 



I. INTRODUCTION 

The quest for novel states of matter remains a central 
theme in physics and one that spans all of its subfields. 
The discovery of high-temperature superconductors and 
the experimental realization of Bose-Einstein condensa- 
tion are some recent examples. The widespread impact of 
these discoveries in the areas of condensed- matter physics 
and atomic physics is well documented in the literature. 
At the interface between nuclear and particle physics — 
but with impact in cosmology and astrophysics — is the 
search for the quark-gluon plasma (QGP), a novel state 
of matter that may have existed at the dawn of the Uni- 
verse and that may still exist today in the dense envi- 
ronments of exotic stars. The quark-gluon plasma, a de- 
confined state of quarks and gluons, is predicted from 
QCD to be attained at high temperatures and/or high 
baryon densities. A variety of experiments have been 
devoted to produce the QGP in the laboratory. These 
efforts started at the AGS in Brookhaven and at the SPS 
in CERN, and now continue at the relativistic-heavy-ion 
collider (RHIC), and should culminate with the construc- 
tion of the large-hadron-collidcr (LHC) at CERN Q. By 
colliding extremely energetic heavy ions, the aim of these 
experiments is to create a region of such high-energy den- 
sity (or high temperature) that quarks and gluons will 
become deconfined. While current experimental facilities 
(primarily RHIC) may have already created the QGP [2(, 
a great challenge remains: how to identify clearly and un- 
ambiguously the production of such coveted state? 

Spatial observations complement terrestrial searches 



for the QGP. The advent of sophisticated telescopes op- 
erating at a variety of wavelengths have turned neutron 
stars from theoretical curiosities into powerful diagnostic 
tools. The observation of an anomalous mass-radius re- 
lation and/or an enhanced cooling in neutron stars, may 
provide strong evidence in support of strange stars |3|, 



quark stars 4j, or neutron stars with exotic cores. These 
cores may contain novel states of matter, such as meson 
condensates, strange-quark matter, and/or color super- 
conductors 
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In the case of terrestrial searches for the QGP, there 
is a set of complementary observables which provide in- 
formation at different stages of the experiment. For ex- 
ample, leptons and photons produced at the center of 
the collision arc expected to carry valuable information 
concerning the earliest stages of the collisions, while the 
relative abundance of light-flavor hadrons should reflect 
the freeze-out stage. Evidence for the QGP is also im- 
printed in the response of the medium to the formation 
and propagation of a heavy quark-antiquark pair, where 
Debye screening of the bound state (e.g., J/ip) is ex- 
pected 0, @- In this case it is important to study the 
possible dissolution of the bound state and the evolu- 
tion of the confining potential with baryon and/or energy 
density. While the quest for these experimental signa- 
tures will undoubtedly continue, convincing arguments 
in favor of the discovery of a str ongl y-coupled QGP at 
RHIC have recently been made (lCf. This persuasive 
study, containing a host of valuable references, uses ex- 
isting observables — such as bulk collective flow and jet 
quenching — to justify the claim. 

In the present work we study the efficiency of the nu- 
clear medium in screening the color charges of a heavy 
quark-antiquark pair. That is, we study color screen- 
ing as a function of the baryon density rather than as a 
function of temperature. This distinction is important 
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as our T = formalism is unable to shed light on finite- 
temperature calculations, which themselves appear to be 
in a state of flux. Calculations based on potential mod- 
els predict the J/ip-meson to dissolve at a temperature of 
T~ 1.1 T c [ll| (with T c being the critical temperature for 
the deconfining phase transition). However, recent lat- 
tice results seem to suggest that the JA/j-meson survives 
up to a temperature of T ~ 1.5 T c [T^, Ojj. Moreover, 
the originally strong suppression of the J/ip- meson pre- 
dicted by Matsui and Satz has been recently put into 
question. Indeed, it has been argued that this elusive 
suppression may be offset (at least in part) by statistical 
recombination of J/^-mesons from the initially produced 
cc pairs 0, El E3 • We offer no new insights into 
these important problems. Rather, we consider a simple 
model that dynamically interpolates between a hadron- 
based description at low-density and a quark-based de- 
scription at high-density to study the dynamics of the 
J/V'-meson. In essence, we aim at computing a heavy- 
quark potential as a function of baryon density to com- 
plement (and compare to) lattice studies |18lll9. 20] that, 
while successful at finite-temperature, are unable to sim- 
ulate hadronic systems at finite baryon density. 

To correlate the nuclear-to-quark-matter transition to 
the in-medium properties of the J / ip meson, a string- flip 
model is used |21| • Although we will refer explicitly to the 
J I ip meson henceforth, our approach is applicable to any 
heavy meson (or baryon) added to the hadronic medium. 
The string-flip model belongs to a class of constituent- 
quark models that incorporates many of the phenomeno- 
logical features of quantum chromodynamics (QCD). The 
underlying degrees of freedom are constituent quarks that 
interact via confining strings (or flux tubes) that arrange 
themselves instantaneously to the optimal configuration 

of quarks [H |H |H E3 ■ The strin S" ni P model is thus 
reminiscent of the Born-Oppenheimer approximation of 
wide use in atomic physics where the fast degrees of 
freedom (electrons or gluons) adjust instantaneously to 
changes in the slow degrees of freedom (ions and con- 
stituent quarks). The cornerstone of the string-flip model 
is a many-body potential that: a) is symmetric under the 
exchange of identical quarks, b) confines quarks within 
color-singlet clusters, and c) enables these clusters to sep- 
arate without generating unphysical (long-range) van der 
Waals forces [2(j . We regard these as a minimal charac- 
terization of any realistic quark-based model of hadronic 
matter. To our knowledge, there is no "conventional" 
potential, namely, one that may be written as a sum 
of two-body (or even three-body) forces, that satisfies 
the above requirements. Thus, the string-flip potential 
is truly many body; moving a single quark may affect 
the interaction among all the quarks in the system. The 
many-body potential is obtained by solving a difficult op- 
timization problem [271 ] , as one must decide how to assign 
colored quarks into color-singlet clusters. This "quark- 
assignment" problem is meant to represent the optimal 
configuration of gluonic strings. 

The manuscript has been organized as follows. In 



Sec.llllwe define the string-flip potential and discuss some 
of its most important characteristics. Further, we intro- 
duce a one-parameter variational wave function that, in 
spite of its simplicity, is exact in the low- and high-density 
limits. Indeed, at low density the variational wave func- 
tion favors a strong clustering of quarks into color-singlet 
objects ("baryons"), while in the high-density limit only 
Pauli correlations remain. Finally, medium modifications 
to the properties of a heavy meson are addressed via a 
density-dependent heavy-quark potential. In Sec. 11111 wc 
present our results paying special emphasis to the correla- 
tions between the medium modifications to the properties 
of the J/'i/j-meson and the deconfining phase transition. 
A summary and conclusions are presented in Sec. II VI 



II. FORMALISM 

An exact many-body treatment of QCD at finite 
baryon density remains an unsolved problem. Although 
fresh new ideas continue to emerge, we use a Q CD- 
inspired model to study the screening of color charge by 
the nuclear medium. The main feature of the string- 
flip model employed here is an adiabatic many-body po- 
tential characterized by the arrangement of the gluonic 
strings (or flux-tubes) into a minimum-energy configura- 
tion |2ll I22I ] . Such a many-body potential is symmetric 
under the exchange of identical quarks and is able to con- 
fine quarks within color-singlet clusters ("baryons"). Yet 
the confining force saturates within a baryon, thereby 
allowing the clusters to separate without generating un- 
observed long-range van der Waals forces [2(|. While 
the exact functional form of the potential remains un- 
certain, the simultaneous requirements of quark confine- 
ment and cluster separability are likely to require solving 
some type of quark-assignment problem: which quarks 
belong to which color-singlet cluster? In what follows, 
the most important features of the string-flip model will 
be addressed. A more thorough discussion may be found 
in Ref. H|. 



A. The string-flip potential 

Following a choice introduced in earlier references [2^, 
|2S| . one starts by defining the optimal pairing between 
red and blue quarks as follows: 

A 

V RB = min ^ v (r iR , r jB ) , (1) 

i=l 

where denotes the spatial coordinate of the ith red 
quark and tjb is the coordinate of its blue-quark "part- 
ner" [iB 1 > P(iB) = jB). Note that the minimization 
procedure is over all possible Al permutations of the A 
blue quarks and that the confining potential v is assumed 
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harmonic with a spring constant denoted by k. That is, 
v{T iR ,Y jB ) = ^k{r iR - r jB ) 2 . (2) 

The "blue-green" and "green-red" components of the 
many-quark potential are defined in direct analogy to 
Eq. 

A 

Vbg = mp^t;fr* B ,rj- G ) ' ( 3a ) 

i=l 
A 

Vgr = rrrm^ufr iG ,r., fl ) . (3b) 

i=i 

In this way, the many-body potential of the system is 
obtained by simply adding the three different pairwise 
contributions 

V{r u ...,y n ) = Vrb + V BG + V GR . (4) 

The many-body Hamiltonian describing a system of N 
quarks each with mass m and momentum is now given 
by 

N 2 

H = Y,7&; + V(Ti,...,r N ), (5) 
£ — ' 2m 

?:=i 

Note that the many-body potential defined in Eq. (J3J 
is flavor blind. That is, the optimal pairing among the 
quarks is done according to their color but not their fla- 
vor. In practice, one may imagine collecting all quarks 
into red and blue "buckets" according to their color — but 
irrespective of their flavor — and then searching for their 
optimal pairing. This is then repeated for blue/green and 
green/red quarks. In the course of this work we adopt 
units in which k = TO = 1. 



B. The variational wave function 

The complicated many-body dynamics may be cap- 
tured with a simple one-parameter variational wave func- 
tion of the following form: 

A (r x , ...,r N ) = e ^ AV ( ri - r ")$ FG ( ri , . . . , r N ) . (6) 

Here A is the sole variational parameter, V is the many- 
body potential defined in Eq. and $j?g is a Fermi- 
gas wave function consisting of a product of Slater de- 
terminants. The variational parameter may be regarded 
as the order parameter for the nuclear-to-quark-matter 
transition, as A^ 1 / 2 represents the length-scale for quark 
confinement. Indeed, at low density the average inter- 
quark separation is much larger than the confining scale 
(A -1 / 2 ~ I) and the clustering of three quarks into color- 
singlet clusters (nucleons) ensues. Note that while the 
interactions between quarks within a single nucleon are 
strong, the residual nucleon-nucleon interaction is weak, 



as the color force saturates within each individual nu- 
cleon. This many-body feature of the model precludes 
the development of long-range van der Waals forces. As 
the density increases, the average inter-quark separation 
will become comparable to the confining scale. This will 
signal the transition to the quark Fermi-gas phase. In the 
high-density regime the interactions between quarks are 
weak (asymptotic freedom) and the system "dissolves" 
into a free Fermi gas of quarks. No correlations between 
quarks remain, except those induced by the Pauli exclu- 
sion principle. Further details on the variational wave 
function may be found in Ref. 28] . For reference, the 
value of the variational parameter for an isolated nucleon 
is given by Ao = 1/ \/3 

C. Variational Monte Carlo 

The structure of the variational wave function entails 
an important simplification, as the expectation value 
of the kinetic energy, which involves derivatives, may 
be simplified through an integration by parts. This 
yields [28j ]. 

<*a|T|^ a ) = (T) x = T FG + 2X 2 {W)x ■ (7) 

Here Tfg is the kinetic energy of a (A = 0) free Fermi gas 
and (W}\ is given by 

N 1 2 

w/ = E^( x »-y») - ( g ) 

n=l 

where the sum is over all quarks in the system and y n rep- 
resents the average position of the two quarks connected 
to the nth quark located at x n . The increase in the ki- 
netic energy of the system relative to that of a free Fermi 
gas is the result of clustering correlations. Of course, this 
"penalty" is more than compensated by the potential en- 
ergy which favors the clustering of quarks into nucleons at 
low densities. The expectation value of the total energy 
of the system now reduces to the following expression: 

E(X) =T FG + 2\ 2 (W) x + (V)x . (9) 

This form is particularly convenient because the two 
functions that remain to be evaluated (V and W) are 
local and may therefore be computed using standard 
Monte-Carlo techniques p9j . 

D. The heavy-quark potential 

In this section we introduce the heavy- quark potential 
to correlate the modifications of the in-medium proper- 
ties of the J/^-meson to the deconfinement phase transi- 
tion. The heavy-quark potential is constructed as follows. 
Consider the original system of N quarks described by a 
variational wave function as in Eq. © . We assume that 
the variational parameter A has been fixed at its optimal 
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value. After a large number of thermalization sweeps, 
one measures the expectation value of the potential en- 
ergy [Eq. ©] 



1 M 
lim iVw> 

M^oo M ^ 1 



(10) 



where the M configurations of quarks are distributed ac- 
cording to the square of the variational wave function. 
At the same time that one pauses the Monte Carlo pro- 
cedure to measure the potential energy in the m t h con- 
figuration, one introduces a heavy quark-antiquark (e.g., 
cc) pair separated by a fixed distance r. This procedure 
is implemented in our simulations by adding a triplet of 
static quarks to the system; red and blue quarks fixed 
at the origin and a green quark fixed a distance r away. 
The red and blue quarks are fixed at the same location 
to simulate a source of anti- green color. As the heavy 
quarks are assumed static, their masses play no role in 
the simulations. Having added the three heavy quarks 
into the system, one now evaluates the potential energy 
of the (A + l)-body system by re-computing the optimal 
pairing of quarks into hadrons. (Note that the added red 
and blue heavy quarks are constrained to pair with each 
other.) That is, 



1 M 



,4 m V). (ii) 



The heavy-quark potential is defined as the difference 
between the potential energy of the (A+ l)-body system 
relative to that of the A-body system: 



V QQ (r) = (V) A+ i(r) - (V) A 



(12) 



Note that the heavy-quark potential is computed in the 
sudden approximation. That is, the many-quark wave 
function [Eq. ©] (and thus the location of all the light 
quarks in the system) is assumed to remain unchanged 
as the heavy quarks are introduced into the system. Not 
so, however, the pairing. The presence of the medium 
screens the interaction between the heavy quarks by find- 
ing the optimal pairing of the N — 3(A + 1) quarks into 
hadrons. In this way, the potential energy between the 
QQ pair gets modified relative to its free-space form be- 
cause of the screening. In general, Vnn(r)^kr 2 . 



E. The J/ip meson in free space 

To quantify how a J / tp- meson (or a comparable heavy 
meson) is modified by color screening, we review its free- 
space properties in a model with harmonic confinement. 
Consider a J / ifj-meson as a nonrelativistic system of a 
quark-antiquark (cc) pair of mass m c interacting via a 
harmonic confining potential of spring constant 2k. That 
is, 



H 



2m c 



_P|_ 

2m c 



-(2fc)(n-r 2 ) 2 



(13) 



Introducing center-of-mass and relative coordinates 

R=-(ri+r 2 ), r = (n-r 2 ), (14) 

enables one to reduce the above Hamiltonian to the fol- 
lowing simple form: 



H = 



p2 ( p 2 



4m c \m c 



(15) 



Describing the motion of the system relative to its cen- 
ter of mass, the ground-state wave function of the J/ip- 
meson in the present harmonic approximation becomes 



,~r 2 /2b 2 



*J/V>( r ) = 1XcXc)im , b = (km c ) 



-1/4 



(16) 



Here b is the harmonic-oscillator length and x represents 
a spin-1/2 spinor. Using the above wave function, the 
ground-state energy and mean-square radius of the J /ip- 
meson in free space are given by 



Ej/t}> = 2m 



'J/V 




0.168 . 



(17a) 
(17b) 



The mean-square radius has been obtained by evaluating 
the ground-state expectation value of the operator r 2 /4. 
Moreover, the arrows in Eq. I|17|) indicate the appropriate 
numerical values for the energy and the mean-square ra- 
dius in units in which k = l and m c /m = 5. Note that the 
mass of the charm quark (~ 1.5 GeV) has been fixed at 
five times the value of the mass of the light constituent 
quarks (~ 0.3 GeV). Our goal now is to quantify how 
these baseline values change as the density of the medium 
increases and, further, to correlate these changes to the 
deconfining phase transition. 



III. RESULTS 

Variational Monte-Carlo calculations of a large num- 
ber of ground-state observables, such as the energy, the 
two-body correlation function, and the strangeness-per- 
baryon ratio have been reported in an earlier publica- 
tion j2^. Here we borrow values for the optimized vari- 
ational parameter to compute the heavy-quark potential 
as indicated in Eq. (|12fl • Results obtained from these sim- 
ulations are displayed in Fig. ^ for a variety of densities 
around the deconfining transition density of p c ~ 0.09. 
Recall that the density is given in dimensionless units 
(for a conversion to physical units see Ref. [28]). The 
symbols in Fig. ^ denote the outcome from the Monte- 
Carlo simulations while the lines represent analytic fits 
to the simulation data. While it is evident that quantita- 
tive changes develop below the critical density as a result 
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of screening, no major qualitative changes emerge. For 
example, the QQ potential remains confining for all den- 
sities below the critical density p c . However, at densities 
above p c , a dramatic change in the shape of the poten- 
tial emerges. Indeed, the heavy-quark potential ceases to 
be confining and a QQ spectrum supporting both bound 
and continuum states develops. This behavior is reminis- 
cent of the temperature dependence of the heavy-quark 
potential displayed in Fig. 6 of Ref . 01 ■ 
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mean-square radius (see Fig. yj. For comparison, at a 
density just above the transition (p = 0.10), the energy 
decreases by almost 50% and the mean-square radius in- 
creases by close to 25%, relative to their free-space values. 
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FIG. 2: Ground-state densities for the J/ip meson weighted 
by r 2 /4 for a variety of densities. The quantities in brackets 
are the corresponding mean-square radii, which are obtained 
by integrating the area under the curves. 



FIG. 1: Density dependence of the heavy-quark potential as 
a function of the separation r between the heavy-quark pair. 
The p = potential, Vqq(j-) = r 2 , represents the baseline (free- 
space) potential. 

Medium modifications to the properties of the J/ip me- 
son are studied by solving Schrodinger's equation in the 
presence of the heavy-quark potential. Recall that we 
have fixed the mass of the charm quark (in units of the 
light-quark constituent mass) to m c /m = 5. Expressing 
the "reduced" J/ip wave function as 

*J/*(r) = ^\XcXc)im , (18) 

the Schrodinger equation for the ground-state (I = 0) 
wave function takes the following simple form: 

(-^Js + v qq^ ')) ^ = E r<r) • ( 19 > 

Here fi c = m c /2 is the reduced mass of the cc pair. In 
Fig. we display the ground-state density f 2 (r) weighted 
by r 2 /4, so that the area under the curve yields directly 
the mean-square radius of the J/i/i-meson, which we have 
enclosed in brackets. For densities below the transition, 
the wave functions cluster around the free-space value. 
Indeed, the change in the mean-square radius of the J/ip- 
meson between p = and p = 0.08 amounts to less than 
6%. However, as soon as the transition density is crossed, 
a significant spread in the wave function develops, result- 
ing in both a reduced binding energy and an increased 



Finally, we use Fig. to correlate modifications to the 
properties of the J/ "0-meson in the medium to the decon- 
fining phase transition. Note that all shown quantities 
are given in units of their free-space values. The density 
dependence of the variational parameter A exhibits an 
abrupt drop at a quark density of about p — 0.09. As 
A~V 2 represents the length-scale for quark confinement, 
this drop signals the onset of the transition to quark mat- 
ter. Recall that for A = the only remaining correlations 
in the many-quark wave function are due to the Pauli- 
exclusion principle [see Eq. ©]. While smaller in mag- 
nitude, the corresponding abrupt changes in the energy 
and mean-square radius of the J/ip are unmistakable. 
This provides confirmation of the long-time paradigm 
that modifications to the properties of the J/tp-meson 
in the medium represents a robust indicator of the de- 
confining phase transition. This in a model that, while 
simple, uses exclusively quark degrees of freedom to de- 
scribe dynamically the deconfining phase transition. 



IV. CONCLUSIONS 

Modification to the properties of the J/ip meson in 
the nuclear medium were computed in a constituent 
quark model of hadronic matter, that confines quarks 
within individual color-neutral clusters, yet allows the 
clusters to separate without generating unphysical long- 
range forces. One of the greatest virtues of this string- flip 
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P 

FIG. 3: Density dependence of the energy and mean-square 
radius of the J/lj) meson. The "jump" in these quantities 
at p ~ 0.09 is correlated to the corresponding jump in the 
variational parameter A at the onset of the nuclear-to-quark 
matter transition. All quantities are measured in units of 
their free-space values. 



quark potential. Operationally, this was effected by 
adding two heavy sources of color charge into the system 
and monitoring how the interaction between the charges 
gets modified by the presence of the medium. Screen- 
ing of the heavy-quark charges by the medium yields an 
effective potential that is drastically different from its 
free-space value above the transition density. Although 
the potential below the transition is also modified by the 
medium, it retains its basic form. In contrast, above the 
transition density the potential ceases to be confining 
and a "heavy quarkonium" spectrum that supports both 
bound and continuum states develops. 

In our model one of the main indicators ("the order 
parameter") of the phase transition is the variational 
parameter of the many-quark wave function. This pa- 
rameter, which represents the length-scale for quark con- 
finement, drops abruptly to (almost) zero at the transi- 
tion density indicating how clustering correlations cease 
to be important as the system evolves into the quark- 
matter phase. This drop correlates strongly to similar 
abrupt changes in the in-medium properties of the J/ip- 
meson. This correlation confirms — in a model that uses 
exclusively quarks and flux tubes (or strings) as the fun- 
damental degree of freedom — the long-standing tenet of 
Matsui and Satz Q of using medium modifications to 
the J I -0-meson as a robust signature for the onset of the 
deconfming phase transition. 



model is the description of the deconfming phase transi- 
tion. Indeed, the evolution of the system from a nuclear 
phase, in which quarks clusters into color-singlet hadrons, 
to a quark Fermi-gas phase is dynamical in the model; 
there is no need for additional ad-hoc parameters to char- 
acterize the transition. 

To quantify how the properties of the J/tp meson get 
modified by color screening we computed the heavy- 
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